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ABSTRACT. This paper proposes an experiment to easily detect radiative heat transfer in the 
microwave range. Following an idea given by Pendry more than a decade ago [1], we show that a 
3D array of tnngsten 2/im radins wires with a 1 cm period makes a low cost material exhibiting 
a snrface plasmon in the microwave range aronnd 2.9 GHz. Snch a heated material shonld 
exhibit an emission peak near the plasmon freqnency well above ambient emission. Analysis 
of the signal detected in the near-field should also be a tool to analyze how homogenization 
theory applies when the distance to the material is of the order of the metamaterial period. It 
could also be give a model to non-local dielectric properties in the same conditions. 

INTRODUCTION 

Since the seminal works of Rytov [2], Polder and Van Hove [3], Agarwal |1], it is known that 
energy density can be enhanced from vacuum energy density when the distance involved is lower 
than the wavelength. This energy density is actually a measure of the density of electromagnetic 
states (LDOS) [5] so that a large increase of the LDOS close to a surface implies an increase 
of the thermal energy density. Such an increase of the density of energy occur when the 
material support surface waves such as plasmon polariton or phonon polariton [HI [Tj. This 
enhancement in the energy density has been recently detected by two types of methods. Above 
a dielectric like SiC, near-field infrared radiation has been detected using near-field radiation 
scattering by a metallic probe leading to a new technique named Thermal Radiation Scanning 
Tunneling Microscopy (TRSTM) [S]. Other experiments measuring a heat transfer between a 
sample and a probe above a metal [9] or a dielectric mm have shown near-field thermal 
radiation enhancement. However, such measurements can hardly be done on the spectral 
domain although some of them seems to have been performed recently mm- 

An alternative to such measurements that necessitate to act at the submicronic scale, is to 
make measurements in the microwave domain. In principle, there are two advantages to work 
at such low frequencies. First of all, if we have for instance a resonance in the Ghz i.e for a 
typical wavelength of 10 cm, this resonance should be seen at distances to the material of the 
order of 1 cm. No needs of specific positioning system is therefore required. Moreover, detection 
system in the GHz are able to detect and follow the electromagnetic field. This means that 
an emission spectrum can easily be achieved. The problem is that resonant surface waves 
usually occur in the visible domain for metals supporting surface plasmon and in the infrared 
domain for dielectric supporting surface polariton. We suggest here to use a 3D centimetric 
cubic array of micrometric tungsten wires. Such an array should produce a surface plasmon 
around a frequency of 3 GHz and exhibit a resonant peak in the density of near-field thermal 



Figure 1: Metamaterial made of a 3-dimensional cubic array of cylindrical wires. The wires 
(here tungsten wires) are separated by a in the three directions and have all the same radius 


electromagnetic energy around this frequency. Note that our work follows the pioneering ideas 
developped at Institut Fresnel in Marseille to study radiative properties of complex shape 
particles. Instead of directly study in the visible or the infrared range scattering by complex 
particles or aggregates, these particles are studied in the microwave range where they size have 
been increased in order to keep the ratio size over wavelength constant [m [151 US]. 

In the following, we remind how to calculate effective optical properties of a cubic array made 
of metallic wires. We propose to build such a structure and then use these optical properties 
in order to calculate energy density above such a structure. We show that a signihcative 
enhancement should be detected above these structures at centimetric distances. 

METAMATERIAL SAMPLE 

Let us consider a structure a cubic array made of metallic wires (Fig]T|). We reproduce here a 
calculation originally made by Pendry [1] but giving more details than in the original article. 
Let us denote the wire radius and a the cubic array period. An elementary wire constituting 
the array and crossed by a current I produces a magnetic held H in the orthoradial direction 
of the wire. The magnetic held, given in any textbook reads 
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Here, B is the magnetic induction, whereas po is the vacuum permeability and r the distance 
to the wire axis. Let us now apply the Maxwell-Faraday equation to an elementary frame 
constituting the array (Fig. [Jj). The circulation of the electric held is equal to the time 
derivative of the hux of the magnetic induction through the elementary frame. Note here that 
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Figure 2: Elementary frame of metamaterial. The frame is a square which side is a. The current 
I circulates in the frame. The magnetic held B created by the current is in the perpendicular 
direction to the frame. The electric held E{a) is the electric held in the right part of the frame 
at a distance a from the left part of the frame. 


the time derivative of a quantity will be obtained by a simple multiplication by —iu since the 
quantity considered will be monochromatic with a time dependence chosen as Let us 

hrst calculate magnetic induction hux <F through an elementary frame 
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So that the induction L of an elementary frame can be written L = fiQa/{27i) In a/r^,. Maxwell- 
Faraday equation around the frame reads 
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where R is the resistance of the wire of length a. In terms of the conductivity, this resistance 
expresses R = aji^'nar^. We now suppose that the wavelength involved is much smaller 
than the array period a so that A -C a, then the electric held in the structure will not vary 
a lot at the scale of the elementary frame. The preceding equation therefore gives us a hrst 
relation between the current and the electric held opening the way to the determination of an 
ehective dielectric constant for the material. Indeed, dielectric constant in a medium relates 
polarization to the electric held and the volumetric current is related to the polarization through 
j = dR/dt = —iuR = —iu:eo{er — 1)E. Integrating this relation on an elementary cube 
constituting the array 

al = —iu:eQ{er — l)Ea^ (4) 

Combining equations ([3]) and (jl]), one obtains the metamaterial ehective dielectric constant 
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The right expression is a Drude expression where the plasma frequency can be identihed with 


n 


2 

P 


2t^(? 

ln(a/r^u) 


( 6 ) 


and the damping frequency with 

R eQa^Vtl 

'y = T = —^ 

Plasma frequency expression show that it is of the order of c/a with a very slow correction 
depending on the ration a/r^. The larger this ratio will be, the smaller will be the plasma 
frequency and the larger will be the plasma wavelength Xp = 27rc/Qp. Resonant phenomena 
occur near the plasma frequency and the plasma wavelength. It is important to have this 
wavelength as large as possible in order that it is larger than the array period. Indeed, near- 
held effects occurs at distances smaller than the wavelength. However, this theory is valid only 
if the medium can be seen as an effective medium which can hardly be the case if the distance 
to the material is of the order of the metamaterial period. 


Such a metamaterial could be achieved with tungsten (W) wires. Very small wires (r^ = 2x 10“® 
m) can be found and sold by companies at a price of around 750 € for a 100 m-long coil. Let 
us consider an array of period a = 1 cm, then flp = 2.58 x 10^° rad s ^ so that Xp = 7.31 cm. 
The losses in this model is obtained from the a Drude model htted on the data measured by 
Ordal et ah mi. 
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where oOp = 9.12 x 10^® rad s ^ and oJt = 8.16 x 10^^ rad s ^ 


UJ{UJ + iUJr) 

Knowing that at low frequency, 
the dielectric constant is related to the electrical conductivity a by = 1 + io'/{ueo), the 

electrical conductivity is then given by cr = eoUp/ur = 5.67 x 10^ m“^. In terms of 7 , 

the losses in the metamaterial are 7 = 8.26x10® rad s“^. As it has been seen in the past and 
as it will be seen in the next section, the LDOS and energy density enhancement is the most 
important near the frequency at which the material support a surface polariton. This happens 
when the dielectric constant approaches -1, that is w = f2p/\/2. At this frequency, the field 
penetration depth 6 will be such that A(n)a;5/c = 1 where n = \/\e) is the optical index. If 
e = —1, n = i so that 6 = a\j2 ln(a/r^))/ 7 r = 2.32 cm. Note that the wavelength where the 
resonance occurs is A = \/2Ap = 10.34 cm that is a frequency v = 2.90 GHz. 


If we want a structure which lateral size is larger than the wavelength and which thickness is 
of a few 5 in order to be considered as an opaque body, we need, as a minimum structure a 
parallelepiped of 20 cm x 20 cm x 5 cm. This structure will contain 2000 elementary cubes 
each of them being constituted of 3 cm of wire (note that each wire is shared between 4 cubes 
and that a cube is constituted by 12 1cm-long wires). This structure will therefore need 60 m 
of wires and can be achieved in a 100 m-long tungsten coil cited above. 


NEAR-FIELD THERMAL ENERGY DENSITY ABOVE THE 

METAMATERIAL 


We now calculate the energy density above a semi-infinite metamaterial made of micrometric 
tungsten wires described in the preceding section. This metamaterial is now described as an 
homogeneous model which dielectric properties obey a Drude model. Metamaterial temperature 
is T and the distance of the metamaterial at which the energy density is calculated is z. We 
used the fluctuationnal electrodynamics theory developped by Rytov [2]. Energy density is then 
given by an integration over plane waves (K, kzv) where K is the parallel wave vector and kzv 
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the wave vector perpendicular to the interface in the vacuum such that + kzv = jc?. 
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where 0(a;,T) = fioj/— 1) is the mean energy of an oscillator at frequency u and 
temperature T. r® and are the Fresnel coefficients given by 
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where kz = \/erk^ — K'^. In expression ([2]), the hrst integral is over wave vector which parallel 
wave vector is lower than ko and corresponds to propagative waves. The second integral over 
parallel wave vector larger than ko describes the contribution of evanescent waves which are 
exponentially decaying from the interface. This last contribution can only be seen in the near 
held at sub wavelength distances. It is important when the imaginary part of the Fresnel 
rehection coefficient is large. This happens in p-polarization near the polariton surface wave 
frequency where eiy = — 1 [7] 

In Fig. ini we have represented the density of energy above the metamaterial for 3 different 
distances. We hrst calculated this quantity at a distance of 10 cm from the semi-inhnite 
metamaterial. Actually, at this distance, the energy density is the same as the one we would 
hnd in far-held. The density of energy increases with the frequency around 3 GHz because, 
this part of the spectrum is situated in the Rayleigh-Jeans approximation of the Planck’s law 
where the density of energy increases as the square of the frequency. 

At 1 cm from the interface, the density of energy is completely diherent from the preceding 
case, u exhibits a peak around 3 GHz. We see that around the peak, the density of energy 
has increased by more than one order of magnitude. The peak is situated at 2.9 GHz i.e. for 
u = 1.82 X 10^° rad s“^. It has a width of the order of the metamaterial damping parameter 
7 . The resonant frequency correspond to the frequency where the reflection coefficient in p- 
polarization diverges that is when the dielectric constant approaches - 1 . At this frequency, the 
metamaterial exhibit a surface polariton which is here a surface plasmon. In the same way it 
has been seen for materials supporting surface polariton such as polar material like SiG or Si02, 
the metamaterial exhibits a peak in the density of energy that is also a peak in the density of 
electromagnetic states. 

If the density of energy is measured at a distance of 1 mm, we note that the peak is greatly 
enhanced and that the density of energy has increased by three orders of magnitude. At this 
distance the theory used here is probably note valid because the distance to the material is 
much smaller than the distance between the wires constituting the metamaterial. At this 
scale, one can hardly see the metamaterial as a continuous medium. Indeed, at 2.9 GHz, the 
wavelength involved is of the order of 10 cm which can be considered as large compare to the 
material period. This means that when the plane waves involved in the calculation are mainly 
propagative waves (i.e waves for which K < u/c), spatial variation along the interface are of the 
order oi 2 tt/K > A and larger than the wavelength. When one approaches the metamaterial, 
evanescent waves begin to enter into play. Typically, the largest parallel wave vector involved 


5 








are of the order of 211 jz. This means that spatial variation along the interface will also be of 
the order of 2 ;. Therefore, ^ ~ a is the smallest distance for which one can consider that all the 
contributing evanescent waves have a lateral spatial variation smaller than the metamaterial 
period. Therefore, one can consider our calculation as a plausible approximation as long as 
z > a. When z is smaller than a, it is hard to say at which distance the effective medium 
theory will lose its pertinence. Note however that as long as the theory applies, i.e. for 
distance as small as 1 cm, the effective medium calculation is robust and is not very sensible 
to the metamaterial irregularities. The reasoning leading to the metamaterial optical response 
presented above can be written in any elementary frame even if the material is not strictly 
periodic. What is important is the average wire radius as well as the average frame size a. 
A variation of a few percent in the periodicity a should affect the optical response parameters 
in the same orders since the variations with a and follow a power or logarithmic law. 

We think an experiment detecting the energy density at 3 GHz could help to determine in 
which conditions effective medium theory could apply. This experiment would consist in a small 
antenna or detector that would be approached to the metamaterial and that would measure 
the square modulus of the electric held that is a quantity proportional to the electromagnetic 
density of energy. We also think this experiment could help to understand what are the limits 
of huctuationnal electrodynamics with a local approximation for the materials modeling at very 
short distance. Indeed, for local media, i.e for media which do not exhibit spatial dispersion, 
the energy energy density diverges as 1/z^ when one approaches a heated body. Moreover, 
when two semi-inhnite heated bodies are put at small distance one from each other, it can be 
seen that the 2 ;-dependence in the heat transfer between the bodies is 1/z"^. This divergence is 
of course unphysical for extremely short distance where non-local effects come into play. Such 
local effects have been studied in the case of metals [TH]. The Thomas-Fermi length acts as 
a cut-off length preventing the divergence at subnanometric distances. For dielectrics, it has 
been suggested that the lattice spacing plays a similar role [1^ . The experiment measuring the 
increasing of the electromagnetic energy density should show at which distance the increasing 
stops and where the saturation begins. This should give informations about the scale below 
which local approximation cease to be valid. One could then extrapolate this behavior to the 
case of usual matter and see if it is relevant to set the locality length to lattice constant. 

CONCLUSION 

We have shown in this article that a metamaterial consisting of a regular centimetric cubic 
array of 2 micrometer-radius tungsten wires exhibits surface polaritons in the microwave range 
around 2.9 GHz. Such a metamaterial, that can be build at a low cost price, will exhibit 
when heated a peak in the near-held thermally emitted density of enegy at a typical 1 cm 
distance from the material. These near-held experiment are in principle much more easier to 
realize than in the visible rage because they do not need sophisticated nanopositionning system. 
Making a study of the density of energy variation with the distance should help to understand 
two questions that remain without answer at present time : at what scale the homogeneous 
approximation to describe a metamaterial remains valid? At which distance from a material 
local approximations describing bulk materials are still valid and what is the local length that 
enter into play at extreme distances from a material. 
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Figure 3: Density of energy versus frequency above a metamaterial made of 2 micrometer radius 
tungsten wires with a 1 cm period. The density of energy is calculated for three distances, 10 
cm (top hgure), 1 cm (middle hgure) and 1 mm (bottom hgure). At small distance a peak 
appears in the density of energy at the frequency where the metamaterial exhibit a surface 
polarities (here 2.9 GHz). 
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